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A THEORY  OF  TURBULENCE  IN  A HOMOGENEOUS  FLUID 
INDUCED  BY  AN  OSCILLATING  GRID 

By 
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Several  investigators  have  run  experiments  to  study  turbulence  in  a con- 
tainer of  homogeneo  us  fluid  caused  by  a grid  extending  horizontally  across  the  box 
and  oscillating  vertically.  The  more  recent  experiments  indicate  that  the  root-mean- 
square  horizontal  velocity  <ru  is  proportional  to  the  oscillation  frequency  f and  is 
inversely  proportional  to  distance  from  the  grid,  and  that  the  integral  length  scale  Si 
is  proportional  to  this  distance  and  independent  of  f.  A theoretical  model  is  based 
on  the  Navier-Stokes  equations  and  the  no- slip  boundary  condition  but  contains  no 
other  assumptions.  The  statistically  steady  behavior  is  as  observed  in  the  experi- 
ment. In  the  unsteady  case,  the  front  separating  turbulent  and  non-turbulent  fluid 
propagates  at  a speed  inversely  proportional  to  the  square  root  of  time. 

The  problem  has  importance  with  respect  to  the  mixed  layer  in  the  upper 
ocean  or  th  e lower  atmosphere.  The  present  model  does  not  consider  buoyancy 
effects  but  If  a passive  contaminant  is  present  the  mean  concentration  is  governed 
by  the  one-dimensional  diffusion  equation  with  constant  turbulent  coefficient  of 
diffusion. 


Bouvard  and  Dumas  (1967)  constructed  an  experiment  with  configuration 
represented  in  figure  1.  A grid  is  oscillated  vertically  with  stroke  S and  fre- 
quency f in  a container  of  homogeneous  fluid.  The  motion  induced  is  turbulent 
and  interest  is  attached  to  the  root-mean-square  horizontal  velocity  ctu  and 


the  integral  length  scale  X as  functions  of  distance  from  the  grid  (or  of  distance  z 
from  some  virtual  energy  source),  the  frequency  f,  the  stroke  S and  the  lengths 
Mj, , M3, . . . characteristic  of  the  grid  geometry  and  position.  The  problem  was 
also  investigated  experimentally  by  Thompson  and  Turner  (1975)  and  very  recently 
by  Hopfinger  and  Toly  (1976).  Thompson  and  Turner's  data  suggest  that  ou£  is 
independent  of  distance  from  the  grid  and  that  ctu  is  proportional  to  the  frequency 
of  oscillation  of  the  grid  f at  any  given  level.  These  two  behaviors  were  also 
found  by  Hopfinger  and  Toly.  They  also  measured  i and  found  z,  so  that 
<JU  “ z-1 . 

To  obtain  an  understanding  of  this  experiment,  we  solve  the  problem  of 
statistically  steady,  finite  motions  in  a Newtonian  fluid  filling  infinite  space 
induced  by  a steady  disturbance  source  on  the  plane  z = 0 whose  properties  are 
independent  of  x and  y.  The  properties  of  the  motion  should  be  the  same  as  those 
in  the  experiment  if  the  frequency  of  oscillation  of  the  grid  is  large  and  the  stroke 
and  its  geometrical  lengths  are  small.  A complete  comparison,  however,  requires 
relations  between  some  parameter  or  parameters  characterizing  the  theoretical 
source  and  f,  S,  v,  , M3 

We  begin  by  considering  a thin  "source"  layer  -d  < z < d.  The  layer 
may  contain  a grid  or  some  other  arrangement  conveniently  constructed  in  the 
laboratory  or  in  the  mind.  Integrating  the  averaged  energy  equation  from  d to 
we  get 

0 = -F3  + Es  - Da 

where  F3  and  E3  are  the  energy  fluxes  w(p/p+  qa/2)  averaged  over  the  planes 
distances  z and  d respectively  from  the  plane  z = 0,  and  D3  is  the  total  energy 
dissipation, 


z, 

(1) 

at 
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D^  = -v  I*  (utf3  u+  v^  v + wvaw)dz  (2) 

d 

In  the  expressions  above,  u,  v,  w are  the  velocities  in  the  x,y,  z directions,  p is 
pressure,  p is  the  uniform  density  and  q is  the  fluid  speed. 

If  d is  finite  and  non-zero,  E3  is  finite  and  non-zero  if  we  require  that 
the  velocities  at  all  z be  finite  and  non-zero.  Then,  for  example,  cru  will  vary 
in  some  complicated  way  with  E3,  d,  z,  v and  other  quantities  that  depend  on  the 
nature  of  the  energy  source.  Evidently  we  may  improve  the  situation  by  letting 
d - 0.  Let  us  assume,  tentatively,  that  we  have  the  freedom  to  construct  an 
energy  source  that  produces  disturbances  near  z = d with  the  following  properties: 

(a)  The  only  length  scale  is  d,  i.  e.  derivatives  with  respect  to 

x,y,  z vary  as  d-1  when  operating  on  any  non-constant  quantity. 

(b)  The  velocity  varies  as  cT",  where  n is  a universal  constant. 

(c)  The  vorticity  £ is  proportional  to  the  velocity  scale  divided 

by  the  length  scale,  i.  e. , dwrwl . 

(d)  The  average  energy  flux  on  the  plane  z = d is  of  order  d-3\ 

We  will  see  later  that  we  have  this  freedom. 

Since  Fa  is  finite  and  non-zero,  D2  cannot  be  zero  as  d - 0.  If  D3  is 
finite  and  non-zero,  E3  is  finite  and  non-zero.  But  Ea/d  is  proportional  to  the 
dissipation  function  e just  outside  the  source  layer.  Using  properties  (a)-(d), 
we  get 

— 3 i -3  i 

€ oc  £ oc  d « d 

This  means  that  the  velocities  at  z = d tend  to  zero  so  that  E3  is  zero.  The 
contradiction  reveals  that  E3  and  D3  are  both  large  and  since  F3  is  finite  E3  ~Da. 
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Again,  E3 /d  is  proportional  to  the  dissipation  function  just  outside  the  source 
layer  so  that 

d"30_1  ~ d~an"a 

or  n = 1.  Thus  Ea  tends  to  infinity  but  Ead3  tends  to  a finite  quantity,  which  we 
call  K3.  The  constant  K and  the  viscosity  v are  the  only  two  parameters  of  the 
problems.  Dimensional  analysis  yields,  for  example 

^ =AX  , | =Aa  (3) 

where  Ax  and  As  are  functions  of  the  Reynolds  number  K/v.  As  K/v  - ®,  they 
become  universal  constants. 

We  now  demonstrate  that  a flow  near  z = d may  be  chosen  with  the  four 
properties  listed  above.  We  put  doublets  of  strength  nd3  (finite  h ) at  all  points 
(nd,  md,  0)  where  n and  m are  all  of  the  positive  and  negative  integers.  Doublets 
with  axes  in  the  direction  of  x positive  and  x negative  correspond  to  n+m  even  and 
n+m  odd,  respectively.  We  also  impose  a rigid  plane  at  z = 0 with  infinitesimal 
holes  for  the  doublets.  The  no-slip  condition  insures  that  there  will  be  vorticity. 
The  length  scale  of  the  motion  near  z = d is  d times  a function  of  k/v  and  so  has 
the  property  (a).  The  velocity  scale  near  z = d is  n/d  and  so  has  the  property  (b). 
There  is  vorticity  and  the  scale  is  proportional  to  d“a  so  has  the  property  (c). 
Finally,  an  irrotational  fluid  in  steady  state  has  a uniform  value  of  p/p  + qa/2  and 
it  is  then  clear  that  the  correlation  in  the  mean  energy  flux  term  is  zero.  In  the 
present  problem  there  is  loss  of  energy  because  there  is  vorticity  and  viscosity  so 
that  on  the  average  (p/  p + qa/2)  will  be  less  for  incoming  fluid  than  for  outgoing 
fluid  at  the  plane  z = d.  The  difference  is  Ka/da  times  a function  of  h/v  and  the 


mean  flux  on  z = d is  of  order  h3  /d3,  satisfying  property  (d). 

The  variation  of  ou  and  l with  z in  (3)  agrees  well  with  experiment. 
Equation  (3)  also  reveals  that  cru  is  proportional  to  K.  If  we  now  neglect  the 
effect  of  viscosity,  the  frequency  f is  the  only  experimental  parameter  involving 
time  so  that  cru  should  be  proportional  to  f as  observed.  Notice,  however,  that 
a relationship  between  K and  f is  analogous  to  the  drag  coefficient  in  shearing 
flow  and  the  latter  is  known  to  involve  the  Reynolds  number,  although  weakly, 
when  the  Reynolds  number  is  large  (Monin  & Yaglom,  1971).  Therefore,  <ru/f 
may,  on  close  inspection,  vary  slowly  with  f. 

The  experiments  and  theory  relate  to  the  important  problem  of  the 
mixed  layer  in  the  upper  ocean  and  the  lower  atmosphere,  where,  however, 
buoyancy  is  important  dyna  nically.  The  mean  buoyancy  b,  or  more  generally, 
the  concentration  of  any  contaminant  is  subject  to  the  equation 


where  is  the  eddy  diffusivity.  If  the  buoyancy  variation  is  weak  or,  more 
generally,  if  the  contaminant  is  dynamically  unimportant,  KD  is  proportional  to 
the  constant  eddy  viscosity  and  (4)  may  be  solved  for  b subject  to  initial  and 
boundary  conditions. 

Finally,  consider  the  initial  value  problem  when  the  energy  source  begins 
at  t = 0 in  an  infinite  fluid  at  rest.  A front  separating  turbulent  and  non-turbulent 
fluid  will  move  away  at  a speed  u«,  = © (K,  v,t).  Therefore, 

= A3  (5) 

K* 


where  A3  is  a function  of  K/v  . 
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Figure  1.  Schematic  diagram  of  the  experimental  apparatus. 
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